Graphene-mediated exchange coupling between cobaltocene and magnetic substrates 
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Using first-principles calculations we demonstrate sizable exchange coupling between a magnetic 
molecule and a magnetic substrate via a graphene layer. As a model system we consider cobal- 
tocene (C0CP2) adsorbed on graphene deposited on Ni(lll). We find that the magnetic coupling 
between the molecule and the substrate is antiferromagnetic and varies considerably depending on 
the molecule structure, the adsorption geometry, and the stacking of graphene on Ni(lll). We 
show how this coupling can be tuned by intercalating a magnetic monolayer, e.g. Fe or Co, be- 
tween graphene and Ni(lll). We identify the leading mechanism responsible for the coupling to be 
the spatial and energy matching of the frontier orbitals of CoCp2 and graphene close to the Fermi 
level, and we demonstrate the role of graphene as an electronic decoupling layer, yet allowing spin 
communication between molecule and substrate. 

PACS numbers: 71.15.Mb, 75.50.Xx, 68.43.-h, 81.05.ue 
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The emerging field of organic spintronics capitalizes on 
the novel functionaUties achieved when organic molecules 
are adsorbed on magnetic substrates. The ability to ma- 
nipulate and tune these functionaUties is an important 
goal. Several problems remain however, before such sys- 
tems can be incorporated into new technological devices. 
One in particular is the capability to adsorb molecules 
on surfaces without any detrimental effects being caused 
to either the structural or magnetic properties of the 
molecule. For this reason it is vital to choose molecules 
with maximum structural robustness upon adsorption [1- 
3]. To this end, the phthalocyanine and porphyrin fam- 
ilies are popular choices due to their planar geometry 
[4-9]. However, in some cases, the strong interaction be- 
tween the metal ion of such flat molecules and the sub- 
strate can modify its electronic states and even quench 
its magnetic moment [10]. 

The use of non-planer molecules, such as metallocenes, 
can minimize this effect. Metallocenes are composed 
of a 3d transition-metal ion sandwiched between two 
cyclopentadienyls (Cp). Depending on the metal ion 
species, both non-magnetic and paramagnetic behavior 
can be found [11]. The spin of the metal ion is shielded 
from the surface by the cage formed by the two Cp rings, 
reducing the possibility that it will be modified substan- 
tially after adsorption. Unfortunately, the deposition of 
metallocenes on metallic surfaces is a difficult process [12] 
and, in some cases, complete dissociation of the molecule 
occurs [13, 14]. 

The intercalation of a graphene spacer layer between 
the reactive surface and the metallocene can reduce the 
possibility of molecular dissociation during deposition. 
Additionally, evidence of charge transfer at molccule- 
graphene-Ni(lll) interfaces [15, 16] and the theoreti- 
cal prediction of large charge transfer from cobaltocene 



(C0CP2) to graphene [17] would suggest that a magnetic 
coupling between cobaltocene and the Ni(lll) surface 
through the graphene layer is still achievable. 

In this Letter, we predict, by first principles elec- 
tronic structure methods, a sizable magnetic coupling 
for CoCp2 adsorbed on a graphene layer deposited on 
a Ni(lll) substrate. Furthermore, we propose interca- 
lation of different ferromagnetic metal monolayers, such 
as Fc and Co, between graphene and the Ni substrate 
as a route to tailor the magnetic coupling. Due to the 
unique electronic properties of graphene [18, 19], metal- 
organic systems of this kind could serve as a basis for 
future spintronics devices. 

Density functional theory (DFT) calculations have 
been performed using the projector augmented wave 
method as implemented in the VASP code [20, 21] with 
the Perdew-Burke-Ernzerhof (PBE) exchange correlation 
functional [22]. Dispersion interactions have been in- 
cluded according to the DFT-D2 approach [23]. Further 
computational details can be found in the Supplemental 
Material [24]. For a comprehensive characterization of 
the interface geometry we consider three possible struc- 
tural degrees of freedom, namely the molecular conforma- 
tion, the graphcne-substrate stacking, and the molecular 
adsorption site. Isolated CoCp2 has already been stud- 
ied extensively by DFT [11] and several possible struc- 
tures have been studied. We consider here CoCp2 in the 
D5/1 high-symmetry configuration [24] where two possi- 
ble Jahn- Teller distorted structures characterized by two 
different electronic states occur. The probability den- 
sity of the highest occupied molecular orbital (HOMO) 
of both of these states, labeled ^62 and ^A2, are plotted 
in Fig. 1(a). In both cases, the CoCp2 molecule at- 
tains a nominal S=l/2 spin. The small lattice mismatch 
(1.2%) of graphene and Ni(lll) lattice constant results 




Table I: Total energy difference AE (meV), Co-graphene 
distance d (A), and exchange energy E^x (nieV) for dif- 
ferent structural and electronic configurations of CoCp2 on 
graphene/Ni(lll) for antiparallel alignment of Co and Ni 
magnetic moments. 



Figure 1: (a) Probability density of the CoCp2 HOMO for 
the ■^A2 and ^B2 states, (b) top-fcc and bridge-top stacking 
of graphene on Ni(lll) (the topmost, second, and third Ni 
layers are colored orange, yellow, and grey, respectively), (c) 
Adsorption geometries of the CoCp2 on graphene/Ni(lll) for 
the top-fcc stacking (Co, C, H atoms in CoCp2 are colored 
blue, dark grey, white, respectively). 



in pseudomorphic growth and the flat conformation of 
the graphene layer [25]. DFT calculations have shown 
that the bonding between graphene and the Ni(lll) sur- 
face is primarily due to van der Waals (vdW) interac- 
tions [26], with a binding distance of ^2.1 A, in good 
agreement with experiments [27] . The morphology of the 
graphene/Ni(lll) interface has been investigated experi- 
mentally [28] and two stable configurations observed (see 
Fig. 1(b)). The top-fcc stacking has two inequivalent C 
atoms, one on top of the Ni(lll) surface atom {Ctop), 
the other on the fee site (C/cc)- The carbon atoms of the 
bridge-top stacking are in bridge positions with respect 
to underlying Ni atoms. We find the top-fcc stacking to 
be more stable by 6.5 meV per C atom than the bridge- 
top stacking. Finally, we find the configuration with the 
molecule axis parallel to the graphene layer more sta- 
ble by 40 meV as compared to the case of perpendicular 
orientation, consistent with Ref. 17. 

We have taken into account several possible adsorp- 
tion geometries CoCp2 can assume on graphene/Ni(lll), 
which we labelled as hollow, bridge, and top, depending 
on the position of the Co atom with respect to the C 
atoms below (Fig. 1(c)). The results are presented in 
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Table I, including the total energy difference A£' (with 
respect to the ground state), the Co-graphene distance 
d and the exchange coupling energies -Eea;, defined as 
Eex = Eap — Ep, where Eap {Ep) is the total energy 
of the system when the spin moment of the Co atom is 
antiparallel (parallel) to the one of the Ni slab. Here, a 
negative value of the exchange energy indicates that the 
cobaltocene's spin moment preferentially orients antipar- 
allel to the Ni magnetization. The lowest energy config- 
uration is found when the molecule is adsorbed on the 
hollow site of graphene, which has a top-fcc stacking on 
the underlying Ni(lll) substrate. The calculated adsorp- 
tion energy of this configuration is ^0.64 eV, somewhere 
between those indicating physisorption and chemisorp- 
tion. A comparison of the total energies in Table I 
shows that, except for the case of configuration 2, all 
other configurations are strongly energetically unfavor- 
able. For all configurations, the magnetic ground state 
shows the molecular spin preferentially aligning antipar- 
allel to the Ni magnetisation, with E^x of the order of 
— 10 meV. This energy is remarkably large if we consider 
that the distance between the Co and Ni atoms is approx- 
imately 6.4 A. As a comparison, an exchange energy of 
only 50 meV was found for chemisorbed Fe porphyrin on 
Co(lOO) [4], despite the much smaller Fe - Co distance 
of 3.5 A. Moreover, the values of E^x that we find are 
high enough to ensure the stability of the spin moments 
against temperature-induced fluctuations under typical 
experimental conditions. In all cases, the relaxed adsorp- 
tion distance between the Co ion and graphene lies be- 
tween 4.3 and 4.4 A and, therefore, cannot play a strong 
role in the differing exchange energies. 

To elucidate the physical origin of the molecule- 
substrate exchange coupling we modify independently 
three possible contributions: the CoCp2 electronic state, 
the graphene stacking, and the CoCp2 adsorption site. 
For the first we found that switching from the ^62 to 
the ^A2 electronic configurations (configurations 1 and 
2 in Table I) lowers the exchange energy to —1.3 meV. 
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Figure 2: LDOS of CoCp2 on graphene/M/Ni(lll) in the 
antiparallel configuration, with M = Ni, Co and Fe: (a) 3d 
states of Co of CoCp2 (b)-(d) 3d^2 states of the M layer atoms 
and 2pz states of Ctop and C/cc (for M=Ni) or Ctop and Chcp 
(for M — Co, Fe; see text). The d states of the M layer are 
plotted in grey, while graphene C p states are in orange and 
dark green. Each curve in panels (b)-(d) is the average over 
the three atoms of that species closest to the CoCp2 center. 
Inset in panel (a): 3d states of Co of CoCp2 for the isolated 
molecule. 



This considerable decrease can be attributed to the re- 
duced extent of the CoCp2 spin-polarized HOMO (see 
Fig. 1(a)) which is critical to determining the size of 
the coupling. Varying the graphene stacking from top- 
fcc to bridge-top (configurations 1 and 5) does not influ- 
ence the magnetic coupling in any appreciable way. This 
is somewhat surprising since the magnetic moment in- 
duced on graphene is approximately one order of magni- 
tude smaller in the bridge-top than in the top-fcc stack- 
ing, with values of -t-0.002 /is and -0.03 / +0.02 fiB, 
respectively. We can conclude therefore that the mag- 
netic coupling does not depend on the size of the mag- 
netic moment induced on the graphene atoms. Finally, 
varying the adsorption site (configurations 1,3,4 and 6) 
can change Eg^ by up to a factor of two. However, the 
coupling remains antiferromagnetic in all cases. 

In Fig. 2 we present the local density of states (LDOS) 
of the system in its ground state (configuration 1). Upon 



Table II: Magnetic moments of the two non-equivalent atoms 
of graphene oiq'' ((ib) and m-Q '^^ (/^s)) the interface metal 
monolayer niM (mb), and the exchange energies Eex (meV) 
for CoCp2 on graphene/M/Ni(lll) (M = Ni, Co, Fe). 
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adsorption on the surface, we observe a small shift to 
higher energies of the molecular Co d orbit als with the 
result that the HOMO is pinned to the Fermi level (Ep) 
of the substrate. It also becomes partially depopulated. 
This is accompanied by a charge transfer of 0.28 e~ from 
the molecule to the surface and a decrease of the magnetic 
moment size associated to the Co atom from +0.74 fis to 
+0.47 fj.B- A hybridization between the 2pz orbital of the 
graphene atoms and the 3^22 orbital of the Ni atom is also 
evident, resulting in the polarization of graphene. No- 
tably, only the C/cc atoms exhibits this strong hybridiza- 
tion with the Ni atoms close to Ep. The energy overlap 
between the minority states of graphene and the minority 
d states of CoCp2 just below Ep is responsible for the sta- 
bilisation of the AP alignment. This energy matching is 
absent for the P alignment, due to the inverted HOMO 
spin polarization. We can thus conclude that the spin 
polarisation of graphene close to Ep determines the sign 
of the magnetic coupling. This is further corroborated 
by the analogous situation occurring for configuration 5, 
for which both the graphene LDOS around Ep [24] and 
the magnetic coupling are similar to the ones of config- 
uration 1. Such a dependence suggests that if one can 
modify the induced spin polarization of graphene in this 
energy window, one can modify the magnetic coupling. 
In order to validate this idea, we have performed sev- 
eral additional calculations, intercalating different mag- 
netic monolayers (Fe and Co) between graphene and the 
Ni(lll) substrate. Experimentally the intercalation of Ni 
and Co monolayers between graphene and Ir(lll) [29, 30] 
and of Fe between graphene and Ni(lll) [31] has been 
successfully achieved. As there are no experimental data 
for graphene/Co/Ni(lll), we have used the same struc- 
ture as for graphene/Fe/Ni(lll). As discussed in Ref. 31 
the intercalated Fe atoms are preferentially placed in the 
fee hollow sites of the Ni, following the Ni(lll) stacking. 
On this substrate, graphene adsorbs in a top-hep struc- 
ture, where the two inequivalent graphene C atoms are 
placed alternatively above the Fe atoms and the hep sites 
(corresponding to the topmost Ni layer) [24]. 

We present in Fig. 2 (c) and (d) the LDOS in the 
case of Fe and Co intercalation. The values of the cor- 



responding magnetic moments and the exchange ener- 
gies are hsted in Table II. The magnitude of the spin 
moment in the interface metal (M) layer increases as 
one goes from Ni to Co to Fe and, due to hybridiza- 
tion, this increase also occurs for the moments in- 
duced on the C atoms, i.e., m^^ and rn^'^'^. Counter- 
intuitively, the magnetic coupling is not found to in- 
crease in line with the magnetic moment and in fact 
decreases. We can identify a trend for the exchange 
energy between CoCp2 and the investigated substrate 
from large antiferromagnetic {Ee^ = —9.7 meV) for 
graphene/Ni(lll), to weak antiferromagnetic {Eex — 
—2.3 meV) for graphene/Co/Ni(lll) and weak ferromag- 
netic {Eex = +2.0 meV) for graphcne/Fe/Ni(lll). The 
energy matching between the HOMO of CoCp2 and the 
Pz states of the carbon atoms, which drives the coupling 
between the molecule and substrate, is disrupted by the 
intercalation of the metal layer. The minority d22 states 
of the Co layer lie at higher energies than those of Ni 
and the Fe states are found at even higher energies. Due 
to hybridisation, the Pz orbitals of the graphene atoms 
are similarly shifted to higher energies. This reduces (for 
Co intercalation) and finally prevents (for Fe intercala- 
tion) the energy matching of the C states with the spin- 
polarized HOMO of CoCp2 with a resultant decrease in 
the magnetic coupling. 

Cross sectional plots of the magnetization density, i.e. 
the difference between the spin up and down charge den- 
sities, close to the Fermi level are given in Fig. 3. In 
panel (a) the spatial matching of the CoCp2 HOMO with 
the pz orbitals of the C/cc atom of graphene adsorbed 
on Ni(lll) is evident in the AP alignment where spin 
density lobes from the molecule and the surface atoms 
merge. In contrast, it is absent in the P alignment re- 
sulting in a negative exchange energy. For CoCp2 on 
graphene/Co/Ni (Fig. 3(b)) there is an excess of major- 
ity spin for the Ctop atoms and of minority spin for the 
Cficp atoms which almost cancel each other. However, a 
small preference towards communication through the mi- 
nority spins is suggested by the plot in accordance with 
the weak antiferromagnetic coupling. For the Fe interca- 
lated layer (Fig. 3(c)), the spin density in the graphene 
indicates spin communication for the P alignment, but 
not for the AP, explaining the positive exchange energy. 
The analysis performed for configuration 5 of Tab. I (see 
Fig. 4 in Ref. 24) reveals a scenario similar to the one of 
Fig. 3(a) and is consistent with an exchange coupling of 
similar size. In contrast in configuration 2 spatial over- 
lap between CoCp2 and graphene states at E^? is absent, 
which explains the much weaker magnetic interaction as 
compared to configuration 1. 

Further clarification of the role of graphene in this 
system can be found by considering the situation when 
CoCp2 is adsorbed directly on the Ni(lll) surface. For 
the optimized distance between the Co and Ni atoms 
{d = 4.3 A) we find a charge transfer of 0.64 e~ from 
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Figure 3: Cross-sectional plots of the local magnetization 
density integrated from —0.1 eV to the Fermi level of CoCp2 
on (a) graphene/Ni(lll), (b) gTaphene/Co/Ni(lll), and (c) 
graphene/Fe/Ni(lll). Left (right) panels refer to the parallel 
(antiparallel) configuration. The cross sectional plane, indi- 
cated as purple broken line in Fig. 1(c) (hollow adsorption 
site), cuts through the Co atom and is perpendicular to both 
the substrate and the Cp rings. 



the molecule to the surface which results in a complete 
quenching of the molecular spin. This can be recovered 
by rigidly shifting the molecule away from the surface 
by 1 A whereby the Co ion attains a magnetic moment 
of -1-0.26 hb- As for the case including the graphene 
layer, the magnetic coupling is antiparallel, albeit weakly 
{Eex — —0.4 meV). A further rigid shift of the molecule 
by 1 A results in an increase of the Co magnetic mo- 
ment to -fO.35 fiB while the exchange coupling becomes 
negligibly small. 

In conclusion, our work demonstrates that graphene 
plays a vital role in determining the interaction between 
a magnetic molecule and a ferromagnetic substrate, be- 
having as an electronic decoupling layer, yet allowing spin 
communication. 
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Computational details 



•^^ We performed preliminary calculations on a p(l x 1) hexagonal supercell to obtain geometry optimized 

OQ grapliene/M/Ni(lll) (M = Ni, Co, Fe) structures and a density of states analysis of the systems. In this first 

^_^ part of our work we used a system containing one layer (two C atoms) of graphene and a four-layer metal slab. An 

1—1 optimized F-centered A;-point grid of 17 x 17 x 1 has been used for these small supercell calculations. Graphene layer 

d and the top layer of metal were relaxed along the z-axis. The three bottom layers of Ni were fixed at the bulk geometry 

1-^ with a lattice parameter of 2.49 A that corresponds to the PBE bulk optimized Ni-Ni distance. The geometries of 

(/2 these small cells have been replicated in plane in order to build the bigger p(5 x 5) supercell containing not only the 

1^ graphene-metal slab but also the cobaltocene (CoCp2) molecule. In these second set of calculations, the graphene 

^ and metal layers were fixed and only the CoCp2 was fully relaxed. A F-centered grid of 3 x 3 x 1 k-points has been 

^_i used. Minimization proceeded till forces were lower than 0.01 eV A~^. 

2 We included in our calculations the van der Waals (vdW) dispersion interactions. It is known that the non- 
local exchange-correlation energy functional (vdW-DF) results in an equilibrium distance between graphene and the 
topmost Ni layer larger than 3.5 A [1] for all the calculated structures, in disagreement with experiments. Thus we 

2 employed the semi-empirical potential DFT-D2 of Grimme [2]. With this method the estimate of the interaction 

Q energy is less reliable as compared to the vdW-DF, however the typical errors never exceed 20 % [2]. On the other 

O hand, the DFT-D2 gives trustworthy equilibrium distances. On average the distances obtained with this method 
for small aromatic systems underestimate by about 5 % the experimental values [2] . In light of the fact that in our 

^— I calculations it is essential to get the equilibrium distances as reliable as possible, we decided to use the semi-empirical 

> potential DFT-D2. 

in 

^5 Cobaltocene electronic structure and adsorption on graphene/Ni(lll) 
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Figure 1 




We plot in Fig. 1 the two possible structures of a metallocene molecule, Ds/i and D^d, named after their point 
group symmetry. In the former the two pentagons of the Cp rings are symmetric, while in the latter they are rotated 
by 180°. Total energy calculations show that Dr,h. symmetry is energetically favored for isolated CoCp2 molecules, 
with respect to D^d- The crystal field produced by the Cp rings split the 3d orbitals of Co. If we consider the axis 
of the molecule as our z-axis, the seven 3d electrons are split in: (i) two electrons in a singlet orbital derived from 
d22, (ii) four electrons in a doubly degenerate orbitals derived from dxy and da;2_j,2, (iii) one electron in a doubly 
degenerate orbitals derived from dyz and d^z- The D5/J symmetry is distorted to remove the degeneracy of the frontier 
orbitals, according to the Jahn- Teller effect, so the symmetry is reduced to C2«. There are two possible distortions 
and thus two possible electronic states; depending on whether the Cp rings tilt slightly toward the molecular center 
or outwards, the ^62 or the ^A2 electronic states are respectively produced, as shown in Fig. 1 panel (a) in the main 
text. In the ^62 state, the HOMO has dyz character while in the ■^A2 state it has dxz character. 

The adsorption of CoCp2 on graphene/Ni(lll) can be regarded as strong physisorption as inferred from the fact 
that (i) the smallest distance between the H atoms of CoCp2 and the C atoms of graphene is ~2.4 A, (ii) there is no 
appreciable distortion of the structure of CoCp2 upon deposition on the surface, and (iii) the C - C, C - Ni and C - 
H bonds are non-polar (or only weakly polar) excluding the formation of hydrogen bonds. 

The configuration where the CoCp2 is oriented with its axis perpendicular to the surface and in the top adsorption 
site, i.e. where the Co atom is directly above one of the graphene Ctop atom, is less stable than the configuration 
with the CoCp2 axis parallel to the substrate surface (configuration 1 in Tab. 1 in the main text) by 40 meV . The 
nfinimum distance is 3.0 A, as compared to 2.4 A. 

Graphene stacking for Fe and Co intercalation 




Figure 2 

The surface lattice constant of the Ni(lll) is 2.49 A, with a very small lattice mismatch, i.e. 1.2 %, with the one 
of graphene ( 2.46 A ). In the case of the intercalated Fe or Co monolayer, lattice mismatch increases, yet a single 
layer is stable and does not lead to relevant modification in the graphene structure. The energetically favored top-fcc 
structure for graphene/Ni(lll) (a) compared with the top-hep structure (b) attained for graphene/M/Ni(lll), with 
M = Fe and Co; first (topmost, orange), second (yellow), and third (light gray) metallic layer atoms below graphene 
atoms (dark grey) are depicted in Fig. 2. This figure has been obtained using VESTA [3] 

Analysis of the LDOS and cross sectional plots close to the Fermi level 

We plot in Fig. 3 the spin-polarized LDOS at the C graphene atoms, focussing on the region close to the Fermi level 
(Ep), for configuration 1 and 5 of Tab. 1 in the main text, and the systems with Co and Fe intercalation. The two 
curves in panel (a), (c) and (d) represent the two inequivalent C atoms (Ctop in red, C/cc and Chcp in blue); in the 
configuration of panel (b), only one inequivalent C atom is present. Looking at panels (a) and (b), it is evident that 
close to Ep' spin polarization of the LDOS is similar for top-fcc and bridge-top, although magnetic moments induced 
on C atoms are one order of magnitude apart, i.e. -1-0.020 vs. -1-0.002 fig. This leads to a similar (in size) antiparallel 
magnetic coupling between the molecule and the substrate spins. In panel (c), Co intercalation leads to a reduction 
of the spin polarization; yet the area subtended in the region within 0.1 eV below E^^ in the minority spin channel 
{a —l) is larger than the one in the majority channel {a =t), and the matching with the CoCp2 HOMO results still 
in a (smaller) antiparallel coupling. For Fe intercalation, panel (d), in the same energy region, there is a small excess 
of majority spin, which is consistent with the small parallel coupling between the molecule and the substrate spins. 
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Figure 3 



The Fig. 4 sliows the cross-sectional plot of the local magnetization density (difference between the spin up and 
down charge densities) integrated from —0.1 eV to E^ of CoCp2 on graphene/Ni(lll) for the bridge-top stacking. The 
presence of spin density between the graphene layer and the cobaltocene indicates spin communication for the case 
of antiparallel alignment of the molecule and substrate magnetization. This is missing in the parallel case consistent 
with the negative spin polarization of the LDOS showed in Fig. 3(b). This situation is analogous to the one of the 
top- fee stacking discussed in the main text and explains the similar exchange energy (see Tab. 1 in the main text). 
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